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ABSTRACT 


This  paper  obtains  some  limit  theorems  for  the 
simple  branching  process  allowing  Immigration  (X^)  when 
the  offspring  mean  Is  Infinite.  It  Is  shown  that  there 
exists  a function  U such  that  {e“”u(X|^)>  converges 
almost  surely  and  If  s • xbjlog^U(J)  < -t  where  (bj)  Is 
the  Immigration  distribution*  the  limit  Is  non-defective 
and  non-degenerate  but  Is  Infinite  If  s ■ 

When  s ■ • limit  theorems  are  found  for  {U{X^)) 
which  Involve  a slowly  varying  non-linear  norming. 


Key  words:  B1enaym4-Gal ton-Watson  branching  process. 
Immigration*  martingale  convergence. 


limit  theorems,  regular  variation. 
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1.  IHTRODUCTION  i 

Let  {X^:n«0|l,.. .}  denote  a B1enayae~Ga1 ton-Watson  process  with  | 

Immigration t for  which  the  offspring  probability  generating  function* 
f*  satisfies  m ■ f'(l-)  ■ • and  {I^:n*l*2». . .}*  the  number  of  Immi- 
grants Into  successive  generations,  are  Independent  and  Identically 
distributed  with  probability  generating  function  b satisfying 
b(0)  < 1.  Suppose  also  that  Xq  ■ 0. 

This  paper  describes  the  behavior  of  X^  as  nf«  .In  terns  of  | 

appropriate  limit  theorems.  In  12  we  show  that  for  a certain  Increas- 
ing function  U constructed  as  In  [61.  the  sequence  {e*''u(X^)}  con- 
verges almost  surely  and  that  the  limit  Is  non-defective  If  a certain 
condition  on  the  Immigration  distribution  obtains  and  Is  essentially  | 

Infinite  otherwise.  It  Is  pointed  out  that.  In  the  former  case.  | 

<Xjj)  may  be  classified  In  terms  of  regularity  and  Irregularity  exactly  . 

as  In  the  case  of  no  Immigration  [6]. 

In  13  we  consider  the  case  where  ■ e*”u(X^)  -►  • e.s.  Limit 
theorems  are  obtained  for  which  Involve  a non-linear  norming  by 
a slowly  varying  (SV)  function.  Theorem  2 Involves  no  extra  assump- 
tions but  has  a rather  strange  appearance:  neater  versions  are  given 
In  Corollaries  1-3  under  conditions  expressed  In  terms  of 
6(x)  ■ 1 - b(expI-l/V(e*))) , where  V Is  the  functional  Inverse  of 
U.  These  results  are  analogues  of  the  families  of  limit  theorems  for  | 

the  case  m < ••  which  were  presented  In  Part  I [3].  and  some  examples 
are  given  to  show  that  the  conditions  In  Corollaries  1-3  are  not 

vacuous.  Finally.  It  Is  shown  that  there  Is  no  sequence  (c  ) of  | 

" i 

norning  constants  such  that  Y,/c.  has  a limit  In  distribution  which  i 

n n 

is  neither  defective  nor  degenerate  at  the  origin. 
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To  complete  the  survey  In  Part  I»  It  suffices  to  mention  that  a 

version  of  Theorem  1 Mas  first  proved  In  [2]  under  conditions  on  f 

ensuring  that  U could  be  chosen  as  a suitable  power  of  log^  and  only 
convergence  in  law  was  established.  The  conditions  of  this  result  were 
substantially  relaxed  In  [1]  to  allow  U to  be  SV  at  Infinity  and 
strictly  Increasing.  A restricted  and  more  opaque  version  of  Corol- 
lary 1 was  proved  In  [2]. 

2.  ALMOST  SURE  CONVERGENCE 

Let  f_  denote  the  n^^  functional  Iterate  of  f and  h.  the 
functional  Inverse  of  -log  Let  P„(s)  * n b^exp(-hj|j(s5)| 

and  p(s)  « lim  P (s)  for  0<  s<  -log  q where  q » f (q)  < 1.  Since 

n-^  " 

hj^(s)*0  (n-M.)  It  Is  clear  that  either  p(s)  > 0 or  a 0.  In  the 

former  case  It  Is  clear  that  {Iexp(-hjj(s)X^)]/Pjj(s)}  Is  a martingale 

and  the  martingale  convergence  theorem  shows  that  for  each 

a.s. 

0<s<-log  q,  exp(-h^(s)Xjj)  -*•  W(s),  say,  and  E W(s)  ■ p(s).  The 

behavior  of  {X^}  can  be  analyzed  In  precisely  the  same  way  In  which 
that  of  the  corresponding  process  without  Immigration  was  analyzed  In 
[6].  In  particular  the  classification  of  points  In  (0,-log  q)  as 
regular  and  Irregular  carries  over  to  the  present  case  and  If 
T » sup(s|0  < s < -log  q,  W(s) < 1}  It  follows  that  T Is  non-defective 
and  P(T  zs^)  ■ p{s^),  P(T  =s^)  » 0 for  every  regular  point  s^. 

Finally,  for  the  function  U constructed  In  [6],  the  arguments 
used  In  this  reference  show  that  e“"u(X|^)  *1^*  U(T"^).  In  the  comp- 
lementary case,  when  p(s)  s 0,  the  martingale  convergence  theorem 
shows  that  exp(-h^{s)X^)  0 for  each  0<s<-log  q.  To  proceed 


3. 


further  we  need  the  following  details  about  U.  It  Is  shown  In  [6] 
that  by  starting  with  a fixed  SQc(0,-1og  q)  it  Is  possible  to  con- 
struct U so  that  it  is  continuous  on  tO,-),  identically  zero  on 
lO,l/(-log  q)l,  positive  and  strictly  Increasing  to  Infinity  on 
(l/(-log  q),»)  and  satisfies  the  relations 

(n»0,l,. . . ) 

Choose  a sequence  {Sjjj:in»l,2,. . . } such  that  -log  q>Suj+0  (si-h-). 

Let  be  a subset  of  the  basic  probability  space  such  that 

P(o„)  » 1 and  h„(s„)X„(u))-*-  (n-*-;uen„).  Thus  on  Q*  » H n_» 

M n m n hi  ^ w 

h (s  )X  (u)-^  for  each  m and  this  Implies  that  eventually 
n m n 

X (u)  > l/h.(s„)  and  hence,  from  (1),  eventually 
n n m 

e""u(X^(«))  > U(l/s^)  («en*). 

Letting  n-^  and  them  m-M»,  we  obtain  all  except  the  final  assertion 
of  the  following  theorem. 

Theorem  1.  If  p{s)  >0  then  e""u(X^)  *‘1*  U(l/T)  where  T is. 
defined  above,  and  i f p(s)  ■ 0,  e'’'’u(X^)  •.  Furthcremore « 

p(s)  >0  iff 

I j-1  bj  log^UCj)  < - (2) 

where  (bj)  denotes  the  Immi gratl on  distribution. 

The  last  assertion  follows  from  these  observations.  Since  U 
is  strictly  Increasing,  the  relations  (1)  can  be  solved  to  obtain  an 
explicit  expression 

h„(s)  » l/V(e"u(s‘M)  (3) 

for  h^,  and  the  proof  of  Lemma  2 In  [1]  needs,  with  one  exception. 
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only  trivial  changes  to  show  that  (2)  holds  Iff 

n b[exp(-l/V(e"u(s“^)))]  > 0 

where  V Is  the  functional  Inverse  of  U.  whence  the  assertion  In 
this  case.  The  exception  Is  that  the  proof  given  In  [1]  assumes  that 
log  U Is  SV:  that  this  Is  so  follows  from 

Lemma  1.  If  L Is  SV  at  Infinity  then  so  is  L(U(»)). 

Proof.  By  virtue  of  the  uniform  convergence  theorem  for  SV  functions 
[7]  It  suffices  to  show  that  for  each  x > 0.  U(Xx)/U(x)  Is  bounded 
away  from  zero  and  Infinity  for  all  sufficiently  large  x.  It  suffice* 
to  consider  the  case  1 s x < «.  The  proof  to  be  given  parallels  that 
of  Lemma  2.3.3  In  [6).  Choose  SQc(0»-log  q)  as  in  Construction 
2.3.1  of  161,  and  set  D « [h(sQ),SQ].  The  functions 
9^(5)  * continuous  on  D,  and,  from  equations  (2) 

and  (3)  of  [61,  it  follows  that 

g„{s)  i 9n  + i(s)  ; g„{s)  0 as  n-«-  : (4) 

hence,  by  Oini's  theorem,  this  convergence  Is  uniform  on  D.  Thus  for 
all  n i Oq,  chosen  large  enough,  and  all  ssD,  X/h^(s)  < 

For  each  sufficiently  large  x construct  n and  scO  so  that 
■ hj^(s)  (see  Construction  2.3.1,  [5]).  It  follows  that 

1 <U(xx)/U(x)  < U(l/h^^j(s))/U(l/h^(s))  = e , 

and  the  lemma  now  follows. 

3*  the  case  p(s)  5 0 


The  description  of  the  behavior  of  In  the  case  p(s)  = 0 

Is  through  distributional,  rather  than  almost  sure,  limit  theorems. 


5. 
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The  first  step  In  deriving  them  is  to  note  that,  from  (3), 

E(exp(-h^(t)X^)j  - p^(t)  » ll^bjexp|-l/V(e'"u(t"^)) j J . 

Using  integral  test  comparisons,  as  in  [3],  it  is  not  difficult 
to  show  that 


PpCt)  » A^{t)  exp  log  bj^exp|-l/V(e^U(t"^)) j Jdy 

where  ^ Denoting  the  integral  by 

and  making  the  change  of  variable  e’*  » e^U(t”^),  it  follows  that  if 


t a t -►0+  then 

;n+log  Ud/tp) 

-J-  * I G(x)dx 

dog  Ud/tp) 

v.here  G(x)  • 1 - b|exp (-1/V(e’‘))  j . We  are  now  in  a position  to 
approach  the  main  limit  theorem,  by  making  a suitable  choice  of 
and  by  using  the  following  version  of  Reuter's  Lemma  1 in  [!}• 


(5) 


Lemma  2.  Let  U b^  constructed  accordi ng  to  Construction  2.3.1  in 
[61  and  in  addition  suppose  i t is  stri ctly  increasing  on  (l/(“log  q),«*) 
Let  y^  ^ posi ti ve,  increasing  and  satisfy  the  properties 

(«)  yp(  u)/yj^{  v)-*-»  (n-*^)  i_f  u > v and  0<c<u,  v<di»;  and 

(b)  e”y^(u)-*-  {n-^«;  c<  u<d). 
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Proofs  Let  « {e"”u(W^)  < y„(u)}  • < V(e"yjj(u)) } . Choose 

ue(c.d)  and  c<Uj<u<U2<d.  If  exp|-H^/V(e”y^(u^))|  (i»l,2); 

then  arguing  as  In  [1]  we  obtain 

- exp[-xj^^)  1 P(A^)  < (EyJ^^  exp 

where  x^^^  ■ Y (e”y^(u)) /v (e^y^^Cu^ ))  (i  = l,2).  If  we  can  show  that 

and  0*  assertion  follows  upon  then  letting  -*>u. 

We  Shan  prove  only  that  follows  from  (a)  that  for 

each  t > 1 

jU)  j v(te\(u,))/v(e\(uj)) 

for  all  sufficiently  large  n.  Let  m(n)  be  the  integer  part  of 

n + log  y^Cui)*  By  (**)  o{n)-H»  (n-H»)  and  < 2”yn(u2^  * 

2 

and  hence,  taking  t •*  e . 

xj^^a  V(exp{m{n)+2})/V^exp(m(n)  + l)  (l/V(e))/bjjjjj^j^j(l/V(e))  , 

where  we  have  used  (3).  Hence,  from  (4).  lim  X*^^^  ■ •. 

rlog  X 

Now,  for  X > 1,  let  a(x)  = exp  I G{y)dy,  v/hich  is  SV  at 

■'0 

infinity  and  strictly  Increasing.  Furthermore*  the  ratio 


A(x)/A(xe'’) 


exp 


^log  x>n 

G(y)d: 

nog  X 

exp  (-/®  6(log(zx))dz/z) 


4 


increases  with  x on  (I.-)  from  say,  to  unity,  and  since 

p(t)  = 0 iff  /^G(x)dx  * •,  it  follows  that  (n-^);  moreover, 

for  each  fixed  x > 1,  A(x)/A(xe”)  decreases  to  zero  as  n-*^.  Hence, 
for  any  0<u<l  and  all  n sufficiently  large,  we  may  define 
y (u)>  0 uniquely  in  such  a way  that 
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A(y^(u))  » uA(e"y^{u))  , 

where  U follows  that  Is  Increasing  on  (0,1),  that  y^(u)-^ 

n n 

(n-H«)  and  that  if  0<v<u<l, 


y„(u)  _ ‘‘‘1 

uA(e"yjj(u)) 

A’^{uA(e"y^{v)) 

y„(*)  »■* 

vA(e\(v)) 

A"^|v/.(e"yj^{v)) 

as  n^»,  since  A Is  SV  at  Infinity.  Choosing 

(5)  we  see  that  ■ u and  hence  E|exp|-Xjj/V (e"yjj(u))|  J -►  u (0<u<l), 

and  that  the  other  conditions  of  Lemma  2 are  satisfied.  This  ccapletss 
the  proof  of  the  following  theorem. 

Thereon  2.  If  p(s)=0,  then 

A(e’"u(X^))/A(U(Xj^))  i W 

where  W 1^  unlfo rwly  distributed  on  (0,1 J. 

Analogues  of  Theorem  2 also  obtain  when  n<l  end 
see  [3]. 

The  proof  of  Theorem  2 has  been  carried  out  under  the  condition 
Xq  « 0.  However 

pS"’(t)  • E[e»p(-i!„li„(t))l vi] 

and  since  t was  chosen  to  converge  to  zero,  we  see  that  Theorem  2 
Is  valid  for  any  Initial  state.  Furthermore,  A(e“''u(*))/A(U(*))  Is 
strictly  Increasing  and  continuous  and  hence  the  weakly  convergent 
sequence  In  Theorem  2 is  a Markov  chain  and.  In  addition.  Is  mixing 
[5,  Theorem  2).  Thus  14]  we  cannot  have  convergence  in  probability 
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in  Theorem  2. 

The  quantity  converging  in  law  in  Theorem  2 has  a rather  cdd 
appearance  since  present  in  both  numerator  and  denominator. 

Tidier  versions  can  be  obtained  by  making  estimates  of 


A(e'’^U(X^))/A(U(Xj^))  » exp 


6(zY^)  dz 


(S) 


under  appropriate  assumptions  on  the  behavior  of  6.  For  brevity, 
we  write  = log  U(X^)  throughout. 

Corollary  1.  If  p(s)  = 0 and  lim  xG(x)  = 0,  then 

x-^ 

A(e""u(X^))/A{e")  - H . 
where  W is  uniformly  distributed  on  10,1]. 


Proof.  Since  xG(x)  -►  0,  we  have 

fY 


j " G(y)  dy  = 0 |-log( l-n/Y^)|  ; 
n 


combining  Theorem  2,  (6)  and  (7),  it  fellows  that,  as  n-*-». 


(7) 


n/Y^  -*•  1.  Hence,  as  , 


A(U(X^))/A{e'’)  = exp  I f ''6(y)  dy 


» exp  |o(log[Y^^/n]  )j 


and  the  result  follows. 


Corollary  2.  If  p(s)  = 0 and  limj^^^xG(x)  = a,  0<a<«*,  then 


n"hog  U(X^)  - 1 i W , 
a 


where  PIW  < u]  * {u/(l+u)}  . 


I 

T. 
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Proof.  Immediately,  from  (6),  we  have 

A(e'"u(Xj)/A(U(X^))^<C^'  {1  -n/Y^}*  . 
and  the  result  now  follows  from  Theorem  2. 


Corollary  3,  If  p(s)  = 0,  lim x6(x)  * • , and  G is  regularly 

A ^ * 

varying  at  Infinity,  then  nG(log  U{X^))  ^ W,  where  W is  a standard 
negative  exponential  random  variable. 


Proof. 


Here , 


from  (6) , we  have 
A(e‘"u(X„))/A(U(Xj) 


exp 


-y 

1-n/Y, 


O(zY^)  dz 


;(8) 


combining  (8)  with  Theorem  2 and  lim  ^^xG(x)  = •,  it  follows  that, 
d 

as  n->«»,  n/Y^  0.  Hence,  since  G varies  regularly. 


A(e'"u(X^))/A(U(X^)}  'V  exp  {-nG{Y^)}  . 
and  the  result  follows. 

If  in  Corollary  3 G has  a positive  index  A,  then  0 < A < 1 
and  the  conclusion  can  be  transformed  to  the  form 
a"^  log  U(X^)  ^ W 

where  W has  the  extreme  value  distribution  function  exp(-x“^) 

and  6(a  ) » n"\ 
n 

We  now  show  that  for  any  U the  hypotheses  of  Corollaries  1-3 
can  be  satisfied.  Let  A be  a positive  integer  valued  random  var- 
iable, define  I » rV(A)]  and  let  b be  the  probability  generating 
function  of  I.  Since 

V{A)  -Isis  V(A)  (9) 

and  log  U is  SV  at  infinity,  it  follows  that  condition  (2)  is  satis- 


10. 


fied  iff  E log^A  < «.  In  [3,53.1)  several  examples  were  given  where 
E log^A  = «,  and  for  each  of  these  T(x)  = P(A > x)  is  SV  at  infinity. 
We  now  suppose  this  to  be  always  the  case.  It  follows  then,  from 
Lemma  1 and  (9),  that  P(I  >x)  'v-  T(U(x))  {x->»)  and  hence  an  Abelian 
theorem  for  power  senes  yields 

(l-b(s))/(l-s)  = I s^P{I>j)  -v  (l-s)"^T(U{l-s))  (stl). 
Letting  s = exp (-1/V(e^))  and  invoking  Lemma  1 once  again,  we  obtain 

G(x)  -v  T(e’').  (10) 

Let  logjX  » log  x and  = log(  1 og|^_ jX)  {k=2,3,...)  for 

all  sufficiently  large  x.  In  [3]  an  example  was  given  for  which 

r 

T(x)  'V  d n log.x]  (x-^«) 

k=l 

where  r > 2 and  c is  a certain  constant.  Using  (10)  it  is  ob- 
vious that  /q  G(x)dx  = » and  that  x6(x)  0 (x-^*)  and  hence  this 

example  satisfies  the  conditions  of  Corollary  1. 

Discrete  distributions  were  also  constructed  in  [3]  for  which 
T(x)  a/log  x (0<a<®)  ; 

T(x)  c(log  x)“^'^'^'  (0<c<“,  1<6<2); 

and 

T(x)  'V  (c/b)(log^x)"^  (0<b,c<«,  r > 2). 

Using  (10)  wo  see  that  these  examples  satisfy  the  hypotheses  of 
Corollaries  2 and  3. 

Finally,  we  show  that  the  neatest  version  of  Theorem  2 that 
could  be  hoped  for  is  in  fact  impossible. 

Theorem  3.  There  is  no  sequence  of  constants  (c^)  such  that 
has  a limit  in  distribution  which  is  neither  defective  nor  degenerate 
dt  zero. 


Proof.  For  any  0 < Xj  < Xg  < “.  let  ^ < c“^U(Xjj)  < x^l 

Then,  since  A (e~”u( • )] /A (U( • )]  Is  continuous  and  strictly  increas- 
ing, it  follows  also  that 

Pn(>'l.X2)  = i A(e‘"u{X^))/A(U(X^))  < r^Cx^)!  , 

where 

0 < r^(x)  H ACe'^^c^x)/ a( c^x)  < 1. 

If  (c^)  is  such  that  c”^U(X^)  is  to  converge  in  distribution,  we 

must  have  " because  of  Theorem  1.  Hence,  since  A is  SV, 

we  see  that  [ l+o(  1 ) ] as  n-H».  Choose  any  subsequence 

(n.)  such  that  r (x.)  -►  r for  some  re[0,l].  Then,  for  any  e > 0, 

"k  ^ 

the  intervals  fr„  (x.),r  (x-)]  belong  to  (r-e,r+c)  for  all  k 

i "k  ^ "k  ^ J 

sufficiently  large.  It  follows  from  Theorem  2 that  (x^tX^)  0, 

and  hence  that,  if  c~^U(X^)  converges  in  distribution,  its  limit 
puts  no  mass  on  (0,  ). 

It  is  interesting  to  note  the  contrast  between  the  cases 
p(x)  > 0 and  p(x)  = 0.  In  the  former,  the  asymptotic  behavior  is 
dominated  by  the  underlying  Gal ton-Watson  process,  and  the  effect  of 
immigration,  apart  from  preventing  extinction,  is  seen  only  in  the 
distribution  of  the  limit  of  U(X^)e"'^:  eventually,  the  contribution 
of  the  immigration  process  becomes  negligible.  However,  when  p(x)  =0 
the  immigration  distribution  has  such  a broad  tail  that  U{X^)e”'^  is 
pushed  off  to  infinity  a.s.  by  the  infinite  sequence  of  occasional, 
but  very  large,  inflows  of  immigrants.  The  character  of  Theorem  2 and 
its  Corollaries,  giving  limits  in  distribution  but  not  with  probabi- 
lity one,  reflects  the  nature  of  the  immigration  process  rather  than 
that  of  the  Gal ton-Watson  process.  In  particular,  unlike  the  case 
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when  p(x)  > 0,  the  limiting  distribution  appearing  In  Theoreir.  2 Is 
the  same,  whether  or  not  the  Gal ton-Uatson  process  Is  regular  or 
i rregular. 


h 1i  H It  h 
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